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Abstract 

An operational characterization of the photocycle of the genetic mutant D85N of 
bacteriorhodopsin, BR-D85N, is presented. Steady-state bleach spectra and pump-probe 
absorbance data are obtained with thick hydrated films containing BR-D85N embedded in a 
gelatin host. Simple two- and three-state models are used to analyze the photocycle dynamics 
and extract relevant information such as pure-state absorption spectra, photochemical-transition 
quantum efficiencies, and thermal lifetimes of dominant states appearing in the photocycle, the 
knowledge of which should aid in the analysis of optical recording and retrieval of data in films 
incorporating this photochromic material. The remarkable characteristics of this material and 
their implications from the viewpoint of optical data storage and processing are discussed. 
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Steady-state characterization of bacteriorhodopsin-D85N photocycle 


1 Introduction 

Bacteriorhodopsin (BR) is a chromophoric protein found in the crystalline purple membrane 
(PM) of Halobacterium salinarium, and acts as a light-driven proton pump aiding in the 
organism’s photosynthetic ATP production in oxygen-deficient environments where aerobic 
survival is not possible. For nearly three decades, photochromic films incorporating BR and its 
numerous chemical and genetic variants have been used with considerable success in many 
optical data storage and processing applications; several excellent review articles discuss the 
biophysical and biochemical structure and function of BR as well as its technical applications. 1 
In most, if not all, optical applications, the photochromic nature of BR is exploited to induce a 

suitable absorption and concomitant refractive-index change within the volume of the film, thus 
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accomplishing the desired storage or processing of information carried by optical waves. 

One of the most attractive genetic mutants of BR for optical data storage is D85N, which 
belongs to a class of BR material referred to as blue membrane (BM) and is grown by means of 
the substitution Asp-85 —> Asn in the 248-amino-acid polypeptide chain that forms a 7-segment 

a-helical pocket around the all -trans retinal chromophore. 3 This site-specific mutagenesis 

effectively blocks the proton translocation pathway inside the pocket, hindering the normal 
completion of the molecular photocycle. Thus, unlike the naturally occurring PM (wild-type) 
BR and other popular variants such as D96N, the photocycle of D85N contains an extremely 
long-lived intermediate state that offers the capability of near-permanent data storage. The 
photophysics and photochemistry of D85N, as well as those of the other BM forms of BR that 
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exhibit similar photocycle structures and spectral properties, namely the acidified BM and the 
deionized BM, have been studied extensively; 4-13 our goal here is to provide, based on steady- 
state bleach spectra and pump-probe absorbance data, an operational characterization of the 
BR-D85N photocycle, concentrating on the salient features of this material pertinent to optical 
data storage and processing applications. 


2 Two-state model 

Two-state models have been widely, and successfully, used in the past to explain the optical 
characteristics of BR films. 14 it is therefore tempting to start our consideration of the BR-D85N 
photocycle with the model shown in Fig. 1. Only photo-induced molecular transitions are 
allowed in this model, with possible thermal relaxation of the P state ruled out for the time being, 
as it has been generally assumed to be stable at room temperatures. 

A. Theoretical development 

The photochemical transition rate from state q to state p is designated y l m and can be 
expressed in the form 


r^(A.;?.O=lnlO0„(A o )4^^-/(r,/)=p„(A o )/(r,O. (2-1) 

where (p m ( A) is the quantum efficiency of the q p transition, A q (A) is the absorption 

spectrum of the q state, A 0 is the (vacuum) wavelength of the optical excitation, N 0 is the 
volume number density of BR molecules in the film, d is the film thickness, h and c 0 denote 
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Planck’s constant and the (vacuum) speed of light, respectively, and I(r,t ) is the space- and 
time-varying intensity of the field inside the film. 

Clearly, the absorption spectra of the various molecular states of BR critically affect the 
optical properties of the film. Now, the pure 5-state spectrum A B (X) can be readily obtained via 
a spectrophotometer measurement performed on the unexposed film (in which all BR molecules 
are in the 5 state). Since it is not possible to fully populate the P state, however, the pure 5-state 
spectrum A P (X) cannot be obtained directly, and must be extracted from the measured 
absorption spectrum of the exposed film, which is generally a linear superposition of A g and A P 
with time-varying weights, viz. 


1 ^ 

A(X;t)~ — j[n B (z,t)A B (X) + n p (z,t)A P (X)]dz, ( 2 - 2 ) 

** 0 

where n q (r,t ) denotes the fractional population density of the q state, related to the volume 
number density N q (r,t) of the 4 -state molecules in the obvious fashion as n q = N q j N 0 , and we 

assumed for concreteness and simplicity that the film is bleached by a normally incident plane 
wave, hence resulting in only a z dependence for the intensity and population distributions inside 
the film. Noting that n B + n p =1, the rate equation governing the 5-state population can be 
written down readily as 


dt 


\pl>B (^0 ) + PbP (^0 )] ^ PbP (^ 0 ) A 


(2-3) 
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while, in accord with the Lambert— Beer law, the optical intensity evolves as 

dl lnlO [r ... . „i i 

~ ^ {[A(^o) _ 4>(^o)J rt « + A p (k Q )}l. (2-4) 

Unfortunately, the coupled system of Eqs. (2-3 and 4) cannot be solved generally for n B (z,t) 
and I(z,t ) in closed form, hindering analytical progress with Eq. (2-2). Matters are greatly 
simplified, however, if the film is bleached to steady-state conditions, for which case we readily 
find from Eq. (2-3) as t °c, 0 r equivalently d/dt = 0, that 


n ss o — Pur Oo) 1 

PpAK)+PhA*o)~ i+0(A o )j^’ (2 " 5) 

where the superscript ss 0 denotes steady state reached with bleaching wavelength A 0 , and 
^ = Ppb / Pur ^ the ratio of forward- to backward-transition quantum efficiencies. 

With constant steady-state population densities throughout the depth of the film, implying 
an exponentially decaying steady-state intensity distribution inside the film [cf. Eq. (2-4)], the 

integration in Eq. (2-2) is now trivial, and the steady-state absorption spectrum of the film can be 
written down readily as 


A ss °(A) = n;"A /) (A) + n?A l ,(X) = 


A B a)+<p a 0 ) 


40 . )_r 


(A) 


1 + ^o)^ 


(2-6) 
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where we substituted in Eq (2-5) and made use of n P =1 - n B . Evaluating Eq. (2-6) at A = A 0 
first, we find the hitherto unknown quantity 


A a <P(X 0 )A b (X 0 )A*(X 0 ) 

P ° [l+$a 0 )]d 8 (A 0 )-^a)' 


(2-7) 


We then solve for A P (X) from Eq. (2-6), obtaining, with the help of Eq. (2-7), our final result: 


,( ' [i +<t>(x 0 )}A B a„)-A“-(K) 


We have thus succeeded in expressing the pure P-state spectrum A P entirely in terms of 
the experimentally measured B- and steady-state spectra A B and A". The parameter <&(X 0 ), 
however, remains unknown except for the obvious condition 0 < <Z> < °° imposed by its very 

definition; two different bounds can be derived for this quantity, as described presently. The 
first, and the obvious, requirement is that of positive-semidefiniteness for the calculated pure in- 
state spectrum; i.e., A P (X) > 0 VA, which yields the lower bound 


<*>(A 0 )> 


A~ 0( A 0 ) MD 
AM ^La“»(A). 


- 1 . 


(2-9) 


A second, and perhaps less evident, requirement can now be imposed upon Eq. (2-8) by noting 
that, since the P state can be only partially populated, at wavelengths above the (quasi-)isosbestic 
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point [cf. Fig. 3] the calculated pure P- state spectrum must lie entirely below the measured 
steady-state spectrum. A moment’s reflection reveals that it is necessary to introduce a second 
bleaching wavelength A, > A 0 to be able to bound <P(A 0 ) from above via this argument, and thus 
we have the requirement A P (?i) < A* (A) VA>A„„, which leads to the upper bound 

(/i={A|A£A„}) 


«>a 0 )< 


^a,) inf A„a)-A-\k) 

A,(K) * A^(X)-A‘-W 


- 1 . 


( 2 - 10 ) 


B. Experimental results 

In order to calculate the pure P-state spectrum within the framework of this two-state 
photocycle model, a 2-mm thick BR-D85N film, whose initial (i.e., pure 5-state) absorption 
spectrum is shown in Fig. 2 (initial optical density at 600 nm of iOD 6 oo = 3), 15 was bleached to 

steady state by a collimated beam through a 6-mm x 9-mm aperture, first with a 633-nm He-Ne 

laser, and then with two separate diode lasers operating at 670 nm and 690 nm. During 
exposure, the time evolution of the film transmittance (at the bleaching wavelength) was 
monitored by means of a photodetector to ensure the attainment of steady-state conditions. Once 
in steady state, the film was removed from the bleaching beam and promptly placed in the 
spectrophotometer for measurement of its absorption spectrum. The bleached area was 
irradiated after each exposure with a 442-nm He-Cd laser to restore the film to its initial 
condition before the next experiment. 

Fig. 3 shows the steady-state absorption spectra of the film for these three bleaching 
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wavelengths. The most immediate observation to be made is that larger fractions of BR 
molecules seem to be excited into the P state with bleaching wavelengths that are deeper into the 
red part of the spectrum. T able 1 summarizes the results of the quantitative analysis based on 
these measurements regarding the value of the quantum-efficiency ratio. We note firstly that, 
owing to A 0 > A (as is required for photo-bleaching), the right-hand side of Eq. (2-9) is strictly 
positive, thus providing useful lower bounds on Aq ) for all three bleaching wavelengths. We 

observe also that the further apart the two bleaching wavelengths ( A 0 and A,) are, the tighter the 
upper bound in Eq. (2-10) becomes, asymptotically approaching from above (but never reaching) 
that of Eq. (2-9); indeed, as seen in Table 1 for the case of A 0 = 633 nm, the upper bound on 

0(A O ) obtained for A, = 690 nm is approximately three times smaller than that for A, = 670 nm. 

The distressing observation here, though, is that the lower and upper bounds for any given 
bleaching wavelength A 0 are roughly two orders of magnitude apart, hence failing to provide a 
narrow window for the estimation of <2>(A 0 ). We are thus forced to resort to a somewhat 
arbitrary method of estimation, and we found that, as a rule of thumb, the geometric mean of the 
lower and upper bounds at A 0 provides a fairly satisfactory value for <2>(A 0 ). Since we have two 
distinct upper bounds for Ag = 633 nm, we show the arithmetic mean of the two corresponding 

geometric means in this case, which comes out remarkably close to the experimentally 
determined value ± = 0.0182 quoted in the literature for 630-nm excitation of (aqueous) 

deionized BM. 7 Also, a finite upper bound on <2>(A 0 ) for A 0 = 690 nm cannot be obtained due to 

the lack of a bleaching source with a longer wavelength, and it was found that five times the 
corresponding lower bound gives a good estimate for the quantum-efficiency ratio in this case. 
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As a check for the soundness of our estimation scheme, we find, by linear interpolation based on 
the three values shown in 1 able 1, the value of 0(A O ) at A 0 = 647 nm to be 0.0337, which is 

within the experimentally determined range of = (0.0263, 0.0357) for 647-nm excitation of 

(dried) deionized BM. 3 

We should stress that the quantum-efficiency ratio figures listed in Table 1 are designated 
as typical values due to the ad hoc manner in which they were estimated, and are merely meant 
to exhibit the relative (order-of-magnitude) variation in 0(A O ) as a function of the bleaching 
wavelength. Despite this inevitable lack of precision, it is safe to say that (within the spectral 
region of interest to us) 0(A) seems to increase, apparently linearly and rather dramatically, 
with increasing wavelength. Obviously, this trend may result from an increase in <p PB or a 
decrease in <j> BP (while <j> BP or <p PH , respectively, remain relatively constant), or both taking place 
simultaneously. A more detailed understanding of this behavior requires a complete molecular- 
dynamics simulation at the quantum-chemical level, aided by ultrafast time-resolved absorption 
spectroscopy, to map out the different ground- and excited-state energy surfaces and the 
associated (radiative and nonradiative) transitions in BR— D85N, which is beyond the scope of 
the present effort. 13,16 

With the value of <f>(A 0 ) thus determined, albeit approximately, we turn to evaluate the 
pure P-state spectrum. As A p ( A) represents an inherent characteristic of the material, data for 
different bleaching wavelengths should all yield the same result through Eq. (2-8); indeed, this 
requirement was the primary criterion used above in determining the values of 0(A O ). The 
calculated pure P-state spectra for A 0 = 633 nm, 670 nm, and 690 nm therefore turn out to be 

nearly indiscernible, and Fig. 4 shows the one obtained with 670-nm data, which has the 
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smoothest tail roll-off (i.e., minimal residual spectral contribution from the B state) and was 
therefore judged the most realistic. The steady-state population densities of the two states are 
also shown in Table 1 based on the values of <P(A 0 ) for the three bleaching wavelengths, 
confirming our earlier qualitative observation that longer bleaching wavelengths more fully 
excite the BR molecules from the B state to the P state. 


3 Three-state model 
A. Motivation 

Toward assessing xhe thermal stability of the BR-D85N molecules in P state, a 
pump-probe experiment was performed next whereby the film was first bleached to steady state 
by a pump beam approximately 10 mW in power, 17 and then the film absorbance was monitored 
as a function of time by means of a weak intermittent probe beam, derived from the pump beam 
with a shutter and a 3-OD attenuation provided by a neutral-density filter. (The need for the 
pump and probe beams to be of the same wavelength will become evident in the course of the 
ensuing analysis.) The first 10-minute segment of 633-nm probe absorbance data is given in Fig. 
5, which not only shows that the film absorbance (at the pump-probe wavelength) changes in 
time once the pump beam is turned off (at t = 0), indicating a rearrangement of the steady-state 

population densities through dark (i.e., thermal) processes, but it also reveals that this kinetic 
behavior cannot be modeled adequately with a single exponential (which would presumably 
represent the thermal decay of the P- state molecules), hinting at the presence of one or more 
additional intermediate states in the photocycle. 18 

Further supporting evidence for this is provided by Fig. 6, which shows the evolution of 
the entire absorption spectrum of the film as a function of time after having been bleached to 
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steady state by a 633-nm source, with t = 0 again marking the instant the pump beam is turned 

off. (To facilitate comparison, we mention that the spectra shown in Fig. 3 were all measured 30 
seconds after the pump beams were turned off; i.e., the curve labeled “633 nm” in Fig. 3 and the 
one labeled “30 sec” in Fig 6 are the same.) The key observation here is that although the in- 
state peak at 490 nm appears to be declining very slowly, there is a substantial rise in the 5-state 
peak at 600 nm early on (i.e., within the first 5 minutes). This can be accounted for by the 
thermal decay (back to 5) of a third short-lived (compared to P ) state whose absorbance at 633 
nm is less than that of the B state. 

The same probe-absorbance and time-resolved spectral measurements were then repeated 
with a 670-nm source (data not shown); although a similar trend was observed here as well, the 
increase in film absorbance in the red part of the spectrum and over the same period of time was 
much less prominent in this case. This, in turn, leads to one or both of the following inferences: 
(1) the steady-state population of this third state attained under 670-nm excitation is smaller than 
that for 633-nm excitation, and (2) its absorbance at 670 nm is smaller than its value at 633 nm, 
implying that the peak of its absorption spectrum is likely to be somewhere between those of the 
5 and P states. 

These observations compel us to consider the L state, known to be the dominant 
intermediate state in the BM photocycle, 6 ’ 9-11 as the most likely candidate for the heretofore 
unidentified third state. Smce the P and L states have 9 -cis and 13 -cis retinal configurations, 
respectively, the proposed three-state photocycle model shown in Fig. 7 features a bifurcated 
structure with two competing photocycles 5 <-4 P and 5 <-» L. (The well-known precursor of L, 

the K state, is not included in this model as it is supposed to be very short-lived. 10 ” 12 ) 
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B. Theoretical development 


2 

The population densities of the three states here evolve according to 


n B 7 LB 7 PH 7 HI. + 7 HI. 7bP + 7 BP n H 

n L = 7lh ~7bl ~ 7bi. 0 • n, , (3-1) 

n p y PB 0 ~7bp~7hp _ n p _ 



where y T designates the thermal relaxation rate from state q to state p. Due to the time 
dependence of the photochemical transition rates y p pq [cf. Eq. (2-1)], a closed-form solution 

cannot be written down for Eq. (3-1); fortunately, however, we shall be interested in two separate 
experimental regimes in the sequel, for which Eq. (3-1) takes simpler forms, as detailed below. 

With a sufficiently strong pump beam of wavelength X 0 illuminating the film, one can 
safely assume that y T Bq « y Bq , thus enabling the elimination of the thermal decay rates from Eq. 

(3-1); that this is indeed justifiable will be established subsequently. After adequate exposure, 
the film reaches steady state, and the corresponding molecular population densities can be found 
readily by setting the left-hand side of Eq. (3-1) equal to zero, leading to 







' Ap (A 0 ) 


(3-2) 


where <P q = <p qH / (p Hq are ihe ratios of forward- to backward-transition quantum efficiencies 
between states B and q = L. P. 
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When the pump beam is turned off, which is assumed for convenience to occur at t - 0, 

the molecules in the L and P states are free to relax thermally, thus leading to a redistribution of 
BR molecules among the three states. Since we now have 7 = 0, the photochemical transition 

rates in Eq. (3-1) are identically zero, and thus, using the steady-state solution Eq. (3-2) as the 
initial condition, the population densities are found to rearrange themselves as 


n K (t) 


1 - n,* exp(- f BL t ) - exp(- y T RP t ) 

« (0 

— 

n? exp(-yL') 



n^exp(-r^t) 


for / > 0, where we made use of the fact that n H + n, +n P =\. The film absorbance seen by the 
probe beam at A 0 can thus be written in the form 

A(X 0 ;t) = [A B (X 0 ) A l {1 o) o)]' 

= 4(A 0 )-zi4(A 0 )exp(-7L/)-^4(^o) ex p(-r«pO’ 

/ > 0, where we implicitly defined 


«/.( 0 
«/.(0 


(3-4) 


AA, (A 0 ) = 

AMX 0 ) = n?[MX 0 )-A r (X 0 )\ 
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which give the change in film absorbance due to L- and P-state decay, respectively. Inserting 
from Eq. (3-2) into Eq. (3-5) and solving simultaneously, we obtain 


A B g o pM-AA L (X o)-M>(Ao)K(Ao) 

A 4 ( A 0 )[l + O l ( X 0 )] + ( A 0 )[ A b ( A 0 ) - A A, ( A 0 )] 


and substituting these expressions back into Eq. (3-2), we also find 


AA (A 0 )[H-<i>, .(Aojl+^q,)) A, q„)-4^(A 0 )] 

a,(X,\\+<s>m<,)+*M< o)] 


Finally, as a direct generalization of Eq. (3-4), we express the entire (transient) absorption 
spectrum of the film as a function of time as 


A(X;t) = 4W[l " »r ex p(— - n ? expf-y^t)] 

4 4' l 4(A)exp(-yg / _/) + «^4(A)exp(-y^/). 


(3-8) 


Evaluating this expression at two different time instants t, and t 2 , and then solving the resulting 
system simultaneously for the unknown spectra, we obtain 
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4(A) = 4(*) + 
4(A) = 4(A) + 


^4(A:/ 2 )exp(y SjP / 2 ) 

^4(A;/ 1 )exp(y flP / 

\)-A 


exp (rlpt 

2 )-exp| 

(rlph) 


[exp| 

(rip - rl,. 

>.] 

-exp 

(rl 

-rl } ,]] 

l 


4A;/ 2 )exp| 

( 7 BL ^2 ) 

1 A{X\t x )^x^[y BI 

,)-4 

»w 

exp( 

rlj. 

2 )-exp| 

(rlj^l 

n? \ 

[exp 

(rl,. - rip 

)4 

-exp 

(t'Il 

- rle 

•M 

1 


(3-9) 


Now, the only unknown quantities in the above expressions are the quantum-efficiency 
ratios <P L (X 0 ) and <P p (X 0 ). As was done before, several inequalities such as 


n? > 

4(A 0 )<4(A 0 ) 


4(A) >0 
4(A) >0 


VA 


4(A) < 4(A) 
4(A) < 4(A) 
4(A) < 4(A) 


VA>A r ^ 


(3 -10a) 


(3-106) 


(3- 10c) 


can be imposed upon the calculated quantities in an effort to obtain some useful bounds on these 
parameters; here X red , somewhere between 570 nm and 600 nm, roughly marks the boundary of 
the red part of the spectrum. For the sake of brevity, however, we omit a detailed analysis at this 
point and defer the discussion of relevant results to the next subsection. 

We close here by confirming a point made earlier: if the film is probed for t > 0 at a 

wavelength A, * A 0 , then the absorbance vector in Eq. (3-4) would have to be evaluated at A„ 
essentially introducing two additional unknowns 4(A)) an( * 4(A|) i nt0 the subsequent 


15 



derivation, thus rendering our approach impracticable; hence the need for the pump and probe 
beams to be of the same wavelength, as asserted before. 

C. Experimental results 

Data supporting the three-state hypothesis have already been presented and discussed 
above; what remains to be done is the recovery of the pure /.-state absorption spectrum within 
the present context. Since the changes in probe absorbance seen in 670-nm and 690-nm 
pump— probe experiments are negligibly small (owing to insignificant steady-state L population 
and /.-state absorbance at these wavelengths), reliable parameter extraction and accurate spectral 
recovery can only be done using 633-nm data, with which we shall deal exclusively here as well. 

Eq. (3-4) provides an analytical expression for the absorbance seen by the probe beam as 
a function of time, and thus, by fitting the experimental data with two exponentials [cf. Fig. 5], 
we can readily obtain estimates for the unknown parameters appearing therein, namely the 
thermal lifetimes y T HL and and the absorbance changes AA,(X 0 ) and AA,,(X 0 ) associated 
with the decay of the L and P states, respectively. Before proceeding further, we note from Eq. 
(3-9) that the connection between the pure-state spectra A L (X) and A P (X) and the unknown 
quantum-efficiency ratios <P h (X 0 ) and <P P (X 0 ) is provided via the population density 
expressions in Eq. (3-7). Specifically, we observe that n? is actually a function of 1 + <P P (X 0 ), 
and since it is plausible to expect that & P (X 0 ) = <P(X 0 ) and O(X 0 ) « 1 at 633 nm [cf. Table 1], 
it follows that A L ( X) should depend only very weakly on & P (X 0 ); a similar reasoning holds for 
the dependence of A P ( X) on <P,( A 0 ) as well, and these arguments are indeed found to be true 

when numerical values are substituted into Eq. (3-9). 

This, then, allows us to regard <P P (X 0 ) as known, and attempt to bound <P L (X 0 ) 
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accordingly; in fact, due to our inability to determine it more precisely, we shall simply take 
0 r (X o )= <P(X 0 ) in the sequel, even though there should clearly be a slight difference in the 
value of <f> PB in the two- and three-state models. Eq. (3- 10a) provides the extremely loose 
bounds 4.16 10 -4 < <2>,(A 0 ) < 283.21, whereas Eq. (3-106) yields the considerably tighter 

window 0.0125 < 4 > a (A 0 ) < 20.331. Finally, the third inequality in Eq. (3-10c) provides the 

lower bound <P L (X 0 ) > 0.0160, and we choose to adopt this as a typical value for the quantum- 

efficiency ratio of the 13 -cis cycle, 19 once again our ultimate guide in this admittedly subjective 
selection process being the “visual” quality of the pure-state spectra obtained via Eq. (3-9) using 
these assumed values. 

With all the essential quantities thus at hand, we turn to the evaluation of A,(k) using 
the transient spectra for f, ~ 30 sec and t 2 = 5 min shown in Fig. 6; the result, a somewhat noisy 

spectrum centered roughly around 570 nm, is shown in Fig. 8. As a validation of the 
assumptions elucidated above, we also mention that A p ( A) obtained with Eq. (3-9) is essentially 
identical to the one recovered from the 633-nm steady-state absorption spectrum within the two- 
state formalism. Finally, the values (at 633 nm) of the various quantities pertaining to the three- 
state model are collected ir Table 2, and are to be compared with the corresponding numbers for 
the two-state model given in Table 1 . 

D. The lifetime issue 

A discussion of the thermal lifetimes of the L and P states is finally in order. The initial 
thermal lifetimes T q - ( 7^ ) ' of the two unstable states q = L,P can be estimated on the basis of 
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the theoretical model for the probe absorbance, given analytically in Eq. (3-4) and shown tailored 
to fit experimental data in Fig. 5. This yields the values shown in Table 2, which seem rather 
extraordinary at first glance. 

We note initially that the I-state lifetime in the acidified BM solutions is known to be on 
the order of several milliseconds, 6 ’ 9-1 1 which is roughly three orders of magnitude shorter than 
seen here. One possible suspect for this conspicuous lifetime extension is the gelatin host within 
which the BR-D85N molecules are embedded for the fabrication of the films used in optical data 
storage and processing applications, which may serve to slow down the conformational changes 
in the molecular structure during the relaxation K — » L — » B. 

A perhaps more intriguing observation is related to the R- state lifetime, which is known 
from practical experience to be much longer than found here; in fact, patterns recorded on 
BR-D85N films are retained for months in the laboratory, hence leading to the belief that the P 
state is nearly permanent. 21 This apparent paradox is resolved when we examine the long-term 
633-nm pump— probe behavior shown in Fig. 9, which reveals that the absorbance increase due to 
P- state decay follows a seemingly power-law type trend, leveling off considerably faster in the 
long run than does the exponential time dependence implied by the molecular rate equations. 
(For purposes of comparison, the double-exponential fit shown in Fig. 5 extrapolates to an 
absorbance value of 1.519 after 20 hours, which is considerably higher than the actual value of 
about 1.215 observed in Fig. 9.) 

A convincing hypothesis for this behavior may be based on the fact that the pump-probe 
measurements were performed on a rather large area on the film. It is entirely plausible, then, 
that due to variations in properties like the local pH and the gelatin and water concentrations, P- 
state molecules far enough apart within such an extended volume may exhibit notably different 
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decay dynamics. Consequently, instead of using a constant thermal decay rate y BP for all P-state 
molecules, we now treat this quantity as a random variable / with a suitably chosen probability 

density function p r (y) characterizing the statistics of the fluctuation in the value of y BP among 
the different P-state molecules within the bleached volume of the film. The expected time 
evolution of the P-state population density for t > 0 is then found as 


«,(0 = !«/>('), Pr (7) d 7 

0 

= Jn^ exp(-y/)/? r (y)d y = n s *M r (it ), 

0 


(3-11) 


where M ( co ) = [°° p r (y)cxp(iO)y)d y is the (generally complex-valued) characteristic function 

of y. Although it would be difficult to determine the precise form of p r ( y), some rather general 

requirements can be readily imposed upon it, such as reality, unimodality, causality (i.e., p r (y) 
= 0 for y < 0), and reality and positive semi-definiteness of M r (it). Among the standard 

distribution functions, 21 only the gamma distribution P r (7) = 7fe(f) ^ ' ex p(~f meets 
these criteria, and thus we obtain 


n P (t) = n? 



( 3 - 12 ) 


Here y and p, the mean and the degrees of freedom of the distribution, are parameters at our 
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disposal in obtaining a good fit to the experimental data. We note that with fi -> 0, Eq. (3-12) 

exhibits the desired negative-exponential and power-law behaviors for small and large t, 
respectively, and indeed provides an excellent overall fit for the transient behavior of the film 
absorbance, as shown in Fig. 9 for y = 0.001 1 134 and p = 0.0030431, which is impossible to 

achieve with any pair of negative exponentials. 

Since the P-state population does not seem to decay exponentially, it is not possible, nor 

meaningful, to ascribe to it a ^-type lifetime. In order to quantify the long-term stability of this 
state, we therefore turn to the difference spectra AA(A; f) = A(X;r)~ A(Z;0) shown in Fig. 10, 
measured following the 670-nm bleaching of the film, with A(A; 0) actually taken 30 seconds 
after the pump beam was turned off (i.e., the spectrum labeled “670 nm” in Fig. 3). This figure 
clearly shows the decay of the P- state molecules back to the B state (bleach recovery), and based 
on this data it is estimated that less than 5% of the P-state molecules returned to the B state after 

7 days. 

More interestingly, Fig. 10 also reveals the emergence of yet another state whose 
absorbance is centered around 390 nm, namely the 0 state known to be thermally accessible via 
the P state. 22 Although there is insufficient data to obtain an accurate pure 0- state spectrum, it is 
estimated that nearly 9% of the P-state molecules have transitioned into the 0 state by the end of 

7 days, thus bringing the total loss of P-state molecules within this time period to about 14%. 

With reference to Fig. 10, the P-state molecules were observed to decay back to the B state at an 
initially higher rate, producing only a negligible 0-state population by the end of the first day; 
however, this situation was seen to change with time as 0-state formation started to compete 
with bleach recovery, shifting the isosbestic point between the B and P states from 550 nm to 
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560 nm, with an isosbestic point between the P and Q states emerging around 425 nm. 

This slow formation of Q state is thought to be responsible for the optical irreversibility 
observed with BR-D85N films. 18 Specifically, during each delay (especially during the long 
ones) between spectral measurements after the film was bleached to steady state, a small but 
nonnegligible fraction of the molecules in the P state thermally decay to Q in the dark. The 442- 
nm source used to erase the film after each experiment does not entirely eliminate this residual Q 
formation, and therefore a certain amount of Q - state accumulation is inevitable over a long 
period of time, which translates into a net loss of photoactive material from the photocycle and a 
fatigue of the film. We note parenthetically that a species with a 450-nm absorption peak 
observed in other BM forms of BR was not detected here; 5,13 it is unknown whether this is 
inherent to D85N or due to the presence of the aqueous gelatin host. 

Before closing, we finally tie one last loose end by checking the accuracy of our earlier 
assumption that the thermal decay rates of the L and P states are much smaller than the 
concomitant photochemical transition rates from these states back to the B state, and that they 
can therefore be neglected in analyzing the photo-bleaching of the film. Evaluating Eq. (2-1) 
using the L- and P -state absorption spectra calculated above, reverse-transition quantum 
efficiencies quoted in the literature (<p BI = 0.94, <p BP = 8.8 10“ 3 ), 711 and the optical intensity of 

the 633-nm excitation (7 = 50 w/ m 2 ), we find that y' m is approximately an order of magnitude 

larger than the initial thermal decay rate y T BL corresponding to the ~ e lifetime shown in Table 2, 
and the corresponding rate parameters for the P state are found to be separated by roughly three 
orders of magnitude. It is therefore safe to say that the assumption y' Hq « y p Bq holds indeed. 
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4 Conclusion 


In this paper, we presented a steady-state characterization of the photocycle of the photochromic 
material BR-D85N from the viewpoint of optical data storage and processing applications. 
Based on the simple two-state photocycle model B <— > P and steady-state bleach spectra, we were 

able to recover the pure P- state absorption spectrum, centered at 490 nm, and estimate the ratio 
<P(A) of the forward- and backward-transition quantum efficiencies of this cycle, all in 
remarkably good agreement with experimentally determined results reported elsewhere. It was 
seen in experiments, and also predicted by theory, that longer wavelengths in the red part of the 
spectrum populate the P state more efficiently; it is therefore clear that sources in the deep red 
should be used to obtain the highest possible absorption and index modulations with BR-D85N 
films. Of particular importance is the strong wavelength dependence of which was seen 

to increase roughly linearly with increasing wavelength. This behavior is obviously essential to 
understanding the photochromism of BR— D85N, and seems to be overlooked in several articles 
in the BR literature. 

In trying to assess the thermal stability of the P state, indications for the presence of an 
additional state in the photocycle were seen in pump-probe experiments. A bifurcated three- 
state photocycle model L <-» B P was then proposed and analyzed in light of the pump-probe 

absorbance data and time-resolved spectral measurements obtained from the steady-state- 
bleached film, and we were thus able to reconstruct the pure L-state absorption spectrum, 
centered at 570 nm, within this context. 

The L state seems to decay thermally within a matter of a few minutes entirely back to 
the B state. This observation has a potentially critical implication for data storage and 
processing: clearly, L is formed in regions of high recording-beam intensity where P is in fact 


22 


intended to form; therefore, decay of L back to B after recording is complete essentially 
constitutes a degradation, or contamination, of stored information. Fortunately, as predicted 
theoretically and observed experimentally, the steady-state population density of the L state 
decreases as the excitation wavelength is moved deeper into the red part of the spectrum (e.g., 
less than 0.5% at 690 nm), which also achieves higher P - state populations as mentioned above, 

and therefore this data-contamination issue does not pose a serious problem. 

The P state, on the other hand, was found to be relatively stable over time with a highly 
nonexponential decay behavior. A possible explanation for this was suggested to be the 
randomness in the decay characteristics of P-state molecules within an extended volume, which 
may also be problematic for data storage and processing applications as it implies a non- 
systematic data degradation that cannot be compensated for. If this randomness is in fact due to 
nonuniformities in the pH level and the water and gelatin content of the film, as suggested, then 
this problem can presumably be alleviated with careful film fabrication techniques. 

For very long time delays, difference spectra revealed a thermal product of the P state, 
namely the Q state with an absorption spectrum centered around 390 nm (which could not be 
extracted accurately from available data), and the accumulation of this state over time presents 
yet another dilemma for data storage and processing as it amounts to a net loss of photoactive 
material as well as a degradation of recorded information. It is noted that for repeated 
write-read-erase cycles, it is advisable to irradiate BR-D85N films with a wide-band blue 
source (e.g., filtered white light) to salvage as much of the BR molecules trapped in the Q state 
as possible. For long-term data storage, on the other hand, 0-state formation should be 
monitored periodically, and data should be retrieved or refreshed before the signal-to-noise ratio 
falls below an unacceptable level. 
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Finally, we mention that if this 0-state formation can be enhanced and accelerated by 
means of chemical or genetic engineering, 22 then its presence in the BR-D85N photocycle may 
in fact be turned into an advantage: owing to a further separation of absorption bands, the B « 

Q photochromism would actually provide a larger absorption- and index-modulation depth, 2 
which directly translates into a higher diffraction efficiency (i.e., signal-to-noise ratio) and an 
associated improvement in storage capacity for multiplexed holographic optical memories based 

on BR-D85N. 
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Table Captions 

1 Ratio of the forward- to backward-transition quantum efficiencies and the steady-state 
population densities as a function of the bleaching wavelength in the two-state photocycle 
model of BR- D85N. “Obtained using Eq. (2-9). b Obtained using Eq. (2-10). c See text for 

explanation. ^Obtained using Eq. (2-5) and n r = \ — n H . 

2 Ratio of the forward- 1o backward-transition quantum efficiencies and the steady-state 
population densities for A 0 = 633 nm excitation in the three-state photocycle model of BR- 
D85N. a See text for explanation. b Obtained using Eq. (3-7) and n B = 1 - h, — n P . Based 
on the theoretical fit via Eq. (3-4) shown in Fig. 5. 
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Bleaching 

wavelength 

Quantum-efficiency ratio 

Steady-state 
population densities d 

A 0 (nm) 

A, (nm) 

Lower 

bound" 

Upper 
bound b 

Typical 
value c 

B state 

P state 

633 

670 

0.0018 



0.4854 

0.5146 

690 


670 

690 

0.0045 

0.7878 

0.0594 

0.1845 

0.8155 

690 

— 

0.0491 

oo 

0.2453 

0.0826 

0.9174 


Table 1 
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Quantity 


Typical Value 

Quantum-efficiency 

9 -cis cycle 

®p(K) = 0.0180 

ratios a 

1 3 -cis cycle 

=0.0160 

Steady-state 

B state 

n s * = 0.4733 

population 

L state 

nf =0.0191 

densities b 

P state 

n? = 0.5076 

Initial 

L state 

Tj = 1 min 

thermal lifetimes' 

P state 

Tp = 2.5 days 





Figure Captions 

1 Two-state photocycle model for BR-D85N. 

2 Measured absorption spectrum of the B state. 

3 Mixed absorption spectra of the film in steady state for three bleaching wavelengths; 
isosbestic-point wavelength = 530 nm. 

4 Calculated absorption spectrum of the P state. 

5 Pump ( t < 0)-probe ( t > 0) absorbance data at 633 nm; the theoretical fit employs two 
exponentials corresponding to the thermal decay of the L- and P-state molecules. 

6 Transient absorption spectra of the 633-nm bleached film for three time instants; bleaching 
beam was turned off at t = 0. 

7 Three-state photocycle model for BR-D85N. 

8 Calculated absorption spectrum of the L state. 

9 Long-term pump-probe absorbance data at 633 nm; the theoretical fit incorporates the 
randomness in the thermal decay rate of the P-state molecules. 

10 Difference absorption spectra of the 670-nm bleached film, showing the long-term depletion 
of the P state into the B and Q states. 
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